THE PROBLEM OF DISTINGUISHING BETWEEN A 
CENTER AND A FOCUS FOR NILPOTENT AND 
DEGENERATE ANALYTIC SYSTEMS 



HECTOR GIACOMINli, JAUME GINE^ AND JAUME LLIBRE^ 
We dedicate this paper to the memory of Javier Chavarriga 



Abstract. In this work we study the centers of planar analytic 
vector fields which are limit of linear type centers. It is proved 
that all the nilpotent centers are limit of linear type centers and 
consequently the Poincare-Liapunov method to find linear type 
centers can be also used to find the nilpotent centers. Moreover, 
we show that the degenerate centers which are limit of linear type 
centers are also detectable with the Poincare-Liapunov method. 

1. Introduction and statement of the main results 

Two of the main and oldest problems in the qualitative theory of 
differential systems in is the distinction between a center and a 
focus, called the center problem] and the determination of the first 
integrals in the case of centers, see for instance This paper deals 
with these two problems for the class of analytic differential systems. 

Let p e be a singular point of a differential system in M^. We say 
that p is a center if there is a neighborhood U of p such that all the 
orbits of U\ {p} are periodic, and we say that p is a focus if there is a 
neighborhood U of p such that all the orbits of U\ {p} spiral either in 
forward or in backward time to p. 

Once we have a center at p of a differential system in M^, another 
problem is to know if there exists or not a first integral H defined in 
some neighborhood U of p (i.e. a non-constant function H : U ^ K. 
such that H is constant on the orbits of the differential system), and 
to know the differentiability of H with respect to the differentiability 
of the system. More specifically, we assume that we have an analytic 
differential system having a center at p. Then, it is known that there 
exists a C°° first integral defined in some neighborhood of p, see |25J. 
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It is also known that there exists an analytic first integral defined in 
U \ {p} for some neighborhood U of p, see [21j; but such analytic first 
integral in general cannot be extended to p. For any center p of an 
analytic differential system in it is an open problem to characterize 
when there exists an analytic first integral in a neighborhood of p, or 
simply a local analytic first integral at p. 

A singular point p is a monodromy singular point of a real analytic 
differential system in if there is no characteristic orbit associated to 
it; i.e., there is no orbit tending to the singular point with definite tan- 
gent at this point. Let p be a singular point of an analytic differential 
system. If p is monodromy, then it is either a center or a focus, see 
[T^ ITH] . Moreover, p is a center if and only if there exists a C°° first 
integral defined in some neighborhood of p, see [23] . 

Let p G be a singular point of an analytic differential system 
in M^, and assume that p is a center. Without loss of generality we 
can assume that p is the origin of coordinates (if necessary we do a 
translation of coordinates sending p at the origin). Then, after a linear 
change of variables and a rescaling of the time variable (if necessary), 
the system can be written in one of the following three forms: 

(1) X = -y + Fi{x,y), y = x + F2{x,y); 

(2) X = y + Fi{x,y), y = F2{x,y); 

(3) X = Fi{x,y), y = F2{x,y); 

where Fi{x, y) and F2{x^ y) are real analytic functions without constant 
and linear terms, defined in a neighborhood of the origin. In what 
follows a center of an analytic differential system in R2 is 

called linear 

type, nilpotent or degenerate if after an affine change of variables and 
a rescaling of the time it can be written as system ((T)), Q or Q, 
respectively. 

The characterization of the linear type centers in terms of the exis- 
tence of an analytic first integral is due to Poincare [221 and Liapunov 
[221 5 ^Iso Moussu [2H]- 

Linear Type Center Theorem. The real analytic differential system 
dH) has a center at the origin if and only if there exists a local analytic 
first integral of the form H = x'^+y'^ + F{x,y) defined in a neighborhood 
of the origin, where F starts with terms of order higher than 2. 

An analytic system on the plane will have a singular point of center 
type if a countable number of conditions on the coefficients of the sys- 
tem are satisfied, see 0. Based on the Linear Type Center Theorem 
there is a method, called the Poincare-Liapunov method, which con- 
sists in determining when a system of the form has a local analytic 
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first integral at the origin, and consequently a center at this point. This 
algorithm looks for a formal power series of the form 

oo 

(4) H{x,y) = Y.H^{x,y), 

n=2 

where H2{x, y) = (x^+?/^)/2, and for each n, Hn{x, y) are homogeneous 
polynomials of degree n, so that 

oo 

(5) H = Y.V2u{x^ + y^)\ 

k=2 

where the V2kS are called the Liapunov constants. It is known that the 
Liapunov constants are polynomials in the coefficients of system (^. 
We note that the Poincare-Liapunov method for analytic differential 
systems is an algorithm which at each step uses only a finite jet of the 
system for the calculation of a Liapunov constant. The singular point 
is a center if and only if all the Liapunov constants vanish. For more 
details see and references therein. 

Until now there is no algorithm comparable to the Poincare-Liapunov 
method for determining the center conditions in the case of nilpotent 
and degenerate singular points, except if the singular point has no char- 
acteristic direction because in this last case we can use the algorithm 
of Bautin fTJ (see also P3inil2Z|)- In any case the necessary computa- 
tions for applying Bautin's algorithm are in general more difficult to 
implement that the ones coming from the Poincare-Liapunov method. 
In this paper we shall show that essentially the Poincare-Liapunov al- 
gorithm also works for determining the analytic nilpotent centers and 
a subclass of the analytic degenerate centers. 

Our main result is the following one. 

Theorem 1 (Nilpotent Center Theorem). Suppose that the origin 
of the real analytic differential system ^ is a center, then there exist 
analytic functions Gi and G2 without constants terms, such that the 
system 

(6) X = y + Fi{x,y)+exGi{x,y), y = -ex + F2{x,y) + exG2{x,y), 
has a linear type center at the origin for all e > 0. 

Roughly speaking Theorem Q can be stated saying simply that an 
analytic nilpotent center is always limit of analytic linear type centers. 
Theorem n is proved in Section |21 

By the Linear Type Center Theorem, system (jHI) has a local analytic 
first integral H^{x,y) at the origin for e > 0. If there exists 

\im H,{x,y), 

£\0 
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and it is a function H{x, y) well defined in a neighborhood of the origin, 
then H{x, y) is a local first integral of system ((21) at the origin. Note 
that, in general, H is not analytic, see Remark IT^ 

We note that the Nilpotent Center Theorem reduces the study of the 
nilpotent centers to the case of linear type centers. So, we can apply 
the Poincare-Liapunov method to system looking for analytic first 
integrals of the form H = {ex^ + ?/^)/2 + F{x,y,e), where F starts 
with terms of order higher than 2 in the variables x and y. We deter- 
mine the Liapunov constants from Several examples showing 
the application of the Poincare-Liapunov method to detect nilpotent 
centers are given in Section EJ 

Based in the results obtained for nilpotent centers we establish the 
following definition. 

Suppose that the origin of the real analytic differential system Q is 
a center. We say that it is limit of linear type centers if there exist Gi 
and G2 analytic functions in x, y and e, without constants and linear 
terms in x and y, such that the system 

(7) X = ey+Fi{x,y)+eGi{x,y,e), y = -ex+F2{x,y)+eG2{x,y,e), 

has a linear type center at the origin for all e 7^ sufficiently small. 
A more general definition of limit of linear type centers would be to 
consider functions Gi and G2 that are not analytic in e. 

Theorem 2. Suppose that the origin of the real analytic differential 
system ^ or ^ is monodromy, and that this system is limit of linear 
type centers of the form © 0?" ©, respectively. Suppose also that there 
are no singular point 0/ ® or (|7j) tending to the origin when e tends 
to zero. Then, system ^ or ^ has a center at the origin. 

Theorem |21 is proved in Section |21 The condition that there are 
no singular point tending to the origin when e tends to zero is easily 
verifiable using the lower order terms of the perturbed system (jH)) or 

dZD- 

Another difficulty of the problem of distinguishing between a center 
and a focus becomes from the fact that this problem for degenerate 
centers can be no algebraically solvable:, i.e., it does not exist an infinite 
sequence of independent polynomial expressions involving the coeffi- 
cients of the system, such that their simultaneous vanishing guarantees 
the existence of a center, see El 113 CHI • 

The problem of distinguishing between a center and a focus is alge- 
braically solvable in the class of analytic differential systems of type 
(P) and Q, see Ull 123 1211 • Nilpotent Center Theorem provides 
a new proof of the fact that nilpotent analytic centers are algebraically 
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solvable. This is due to the fact that we have seen in Theorem ^ that 
the nilpotent centers can be e-approximated by systems having hnear 
type centers. Then, applying the Poincare-Liapunov method to these 
hnear type centers, and doing the hmit when e \ we obtain alge- 
braic conditions characterizing the existence of a nilpotent center. So, 
in fact Theorem^ provides an algorithm for solving the center problem 
for nilpotent centers. See the end of Section |21 for more details, and 
also Section m 

For centers of the form Q some preliminary results exist for distin- 
guishing between a center and a focus, see for instance [HI 1111113 US 

We say that an analytic differential system in the plane is time- 
reversible (with respect to an axis of symmetry through the origin) if 
after a rotation 



the system in the new variables {^^rj) becomes invariant by a transfor- 
mation of the form {^,ri.,t) (^, —rj, —t). The phase portrait of this 
new system is symmetric with respect to the straight line ^ = 0. We 
note that for all reversible nilpotent centers which are symmetric with 
respect to a straight line through the origin, this line can be only the 
line of the axes x or y. We remark that all the nilpotent centers that 
we know are time-reversible or have an analytic first integral at the 
origin. 

In the case of degenerate centers is much more difficult to distinguish 
between a center and a focus than in the case of linear and nilpotent 
type centers. In the next theorem we present some results for the 
degenerate centers. 

Theorem 3. For a degenerate analytic center the following statements 
hold. 

(a) A Hamiltonian degenerate center is always limit of linear type 
Hamiltonian centers. 

(b) A time-reversible degenerate center is always limit of linear type 
time-reversible centers. 

(c) There are degenerate centers which are neither Hamiltonian nor 
time-reversible that are limit of linear type centers. 

(d) Non algebraically solvable degenerate centers are not limit of 
linear type centers. 

(e) There are algebraically solvable degenerate centers which are not 
limit of linear type centers. 
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(f) There exist degenerate centers with characteristic directions which 
are limit of degenerate centers without characteristic directions. 

Theorem El is proved in Section 121 

Let Q be a family of analytic systems depending on several param- 
eters. Inside the degenerate centers of this family we can determine 
those which are limit of linear type centers of the form For this 
kind of systems we can apply the Poincare-Liapunov method to sys- 
tem ((7|) with e ^ and compute their Liapunov constants. Vanishing 
these Liapunov constants we obtain the center conditions for the sys- 
tem ((7|). Taking the limit when e — >■ in these conditions, we get the 
center conditions for the degenerate centers System ((7|) must be 
a linear type center only for e sufficiently small. In consequence, 
in the applications of the Poincare-Liapunov algorithm, it is sufficient 
to calculate the Liapunov constants up to first order in e. In contrast, 
for the nilpotent centers, which are always limit of linear type centers, 
we can calculate up to any order in e because system © has a center 
at the origin for all e > 0, and from this fact we can obtain several 
conditions at each step of the algorithm. 

In Section El we provide an example of the application of this method 
to a family of polynomial differential systems Q. 

Finally, in Section El we obtain some results on the cyclicity of nilpo- 
tent and degenerate centers which are limit of linear type centers. In 
particular, we prove the following result (for a definition of cyclicity of 
a center see Sectional). 

Proposition 4. Consider a nilpotent center or a degenerate center of 
a polynomial differential system ^ or ^ of degree m. We suppose 
that this center is limit when e —>■ of linear type centers of polynomial 
differential systems of degree n of the form © or ((7|), respectively. If 
the Liapunov constants of a general perturbation of the same degree n of 
the linear type centers (0) or ((7|) are well-defined when e — > and the 
Poincare map for a perturbation of the initial nilpotent or degenerate 
center is analytic, then the following statements hold: 

(a) The cyclicity of the nilpotent center Q is at most the cyclicity 
of the linear type center © for all e > 0. 

(b) The cyclicity of the degenerate center Q is at most the cyclicity 
of the linear type centers ((7|) for e ^ sufficiently small. 

2. Proof of Theorems [H and (21 

The characterization of the nilpotent centers in terms of the existence 
of a symmetry is due to Berthier and Moussu |^ who obtained the 
following result. We shall need it in the proof of Theorem d 



NILPOTENT AND DEGENERATE CENTERS 



7 



Theorem 5. If the analytic system Q has a center at the origin, then 
there exists an analytic change of variables such that the new system 
has also the form (0) and it is invariant by the change of variables 
{x,y,t) {-x,y, -t). 

We recall from ^Ul that if the analytic system (j2)) has a center at the 
origin and there exists an analytic change of variables such that the 
new system has also the form Q and it is invariant by the change of 
variables {x,y,t) (x, —y, —t), then the system has an analytic first 
integral defined in a neighborhood of the origin. 

Proof of Theorem^ Assume that the origin of system Q is a center. 
Theorem and its proof says that for any nilpotent center Q corre- 
sponding to an analytic vector field X{x,y), there exists an analytic 
change of variables {x, y) — {u, v) of the form 

(8) X = u + . . . , y = V + . . 

such that X{x,y) written in the new variables is a vector field of the 
form 

(9) Y{u,v) = {v + F^{u,v),F2{u,v)), 

where Fi and F2 are analytic functions starting with terms of second 
degree in x and y, and the associated differential system is invariant 
under the change of variables {u,v,t) t-^ {—u,v, —t). 

Now we consider the following perturbation of the vector field Q: 

(10) Y,{u,v) = {v + Fi{u,v),-eu + F2{u,v)), 

with e > 0. Since the eigenvalues at the singular point located at the 
origin are ±y^2, and the differential system associated to the vec- 
tor field ()10|) is invariant under the change of variables {u,v,t) 
{—u,v,—t) (because the unperturbed system is invariant), it follows 
that the origin of the vector field (llUj) is a linear type center for all 
e>0. 

Using the inverse of the change of variables (jH)) we get that the dif- 
ferential system associated to the vector field (|T0|) becomes 

(11) X = y + Fi{x,y)+exGi{x,y), y = -ex + F2{x,y) + exG2{x,y), 

where Gi and G2 are analytic functions without constants terms, de- 
pending on the change of variable (jH)). Let Xe(x,?/) be the vector field 
associated to system (jlip . Since l^e(M, v) has a linear type center at the 
origin for all e > 0, the same holds for Xi,{x,y). This completes the 
proof of the theorem. | 

Proof of Theorem^ Consider an analytic system (P, Q) of the form Q 
or with a monodromy singular point p at the origin. Suppose that 
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this system is limit of linear type centers (P^, Q^) of the form (jUj) or (|7|). 
respectively. Since the origin is monodromy, if is a sufficiently small 
curve with an endpoint at the origin, then the Poincare map U : S ^ S 
associated to the system (P, Q) is well-defined and the leading term 
is always linear for a suitable choice of a semi-transversal algebraic 
curve, which can have a singularity at the singular point, see 1^ - 
The Poincare map Us : S ^ S associated to the system (P^, Q^) is the 
identity for all e > if the center is nilpotent, and for e sufficiently 
small and e ^ if the center is degenerate. Therefore, by the theorem 
on analytic dependence on initial conditions and parameters, it follows 
that n = lim^^o if the center is nilpotent, or 11 = lim^^ong if the 
center is degenerate. Hence, we conclude that 11 is the identity. So, 
the monodromy singular point p of (P, Q) is a center. The condition 
that there are not singular points tending to the origin when e tends 
to zero guarantees that the domain of Us does not reduce to the origin 
when e tends to zero. | 

Theorems n and El can be used to detect nilpotent centers of analytic 
differential systems applying the algorithm of Poincare-Liapunov. In 
the particular case of polynomial systems the method works as follows. 
We consider the system 

(12) x = y + Fi{x,y), y = F2{x,y), 

where Pi and P2 are polynomials without constants and linear terms 
containing a set of arbitrary parameters and such that the origin is a 
monodromy singular point. We recall that using Andreev's Theorem 
we can know when a nilpotent singular point is or not monodromy, see 
[S] . For detecting the centers of ()12p , according with Theorem ^ we 
consider the perturbed system 

(13) X = y + Fi{x,y)+exGi{x,y), y = -ex + F2{x,y) + exG2ix,y), 

where xGi and XG2 are analytic functions starting with quadratic terms 
in X and y. We apply now the Poincare-Liapunov algorithm to deter- 
mine necessary conditions to have a center at the origin for system ()13p. 
In general, these conditions will be satisfied by choosing conveniently 
the coefficients of the analytic functions Gi and G2- When this is not 
possible we must employ the parameters of the polynomial system (jl2|) . 
In this way we will obtain necessary conditions for the existence of a 
center at the origin of system ()12j). The set of sufficient conditions of 
center for the non-perturbed system p2|) will be obtained in a finite 
number of steps, because the Hilbert's basis theorem guarantees that 
this process is finite. Every time that we find a necessary condition for 
the non-perturbed system (fT^ we must to study if the non-perturbed 
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system ()12|1 already have a center at the origin. As the number of steps 
is finite and for determining each Poincare-Liapunov constant of the 
perturbed system (fT^ we need only a finite jet, the necessary pertur- 
bation to detect the center cases will be polynomial, i.e., the functions 
Gi and G2 will be polynomials. We note that, under the assumptions 
of Theorem |2l it is not possible to satisfy the center conditions of (fT^ 
only with the parameters of the perturbation because in that case using 
Theorem 121 the nilpotent polynomial system would have a center for ar- 
bitrary values of the parameters of the family, which is a contradiction 
if the initial system has not a center at the origin. 

3. Proof of Theorem El 

In this section we shall work with an analytic degenerate center Q 
defined in a neighborhood of the origin. 

Proof of Theorem\^a): Suppose that system Q is Hamiltonian with 
Hamiltonian H = H{x,y). The system 

(14) X = -ey + Fi{x,y), y = ex + F2{x,y), 

is also a Hamiltonian system with the Hamiltonian first integral e{x'^ + 
y'^)/2 + H{x, y). Consequently, system (fT^ has a linear type center at 
the origin for e 7^ 0, and the initial degenerate center Q is obtained 
taking in system (jl4p the limit when e — >■ . | 

Proof of Theorem\^b): Without loss of generality, taking into account 
the definition of a time-reversible system, we can assume that system 
Q is invariant by the change of variables {x, y, t) \—>- (x, —y, —t). Con- 
sider the perturbation of it given by a system of the form ()14|) . Then, 
it is easy to see that system (fT^ is also invariant under the change of 
variables (x, t) 1— > (x, —y, —t). Therefore, since the eigenvalues of the 

linear part at the origin of system (fT^ are ±-y/|7[i, it has a linear type 
center at the origin for e 7^ 0. Again, the initial degenerate center Q 
is obtained taking in system ()14|) the limit when e — . | 

Proof of Theorem{^c): Consider the following quartic polynomial dif- 
ferential system 

^^^^ y = (x + 2x2)(x2 + y2). 

It is easy to see that this system has a degenerate center at the origin, 
because removing the common factor x^ + y"^ (doing a change of the 
independent variable) we get a quadratic Hamiltonian system having 
a center at the origin. 



10 



H. GIACOMINI, J. GINE AND J. LLIBRE 



It is easy to check that system (jl5j) is not Hamiltonian, and that it 
has the first integral 

(16) H{x, y) = {x^ + y^)l2 + 2xV3 - y^^- 

Now we claim that system is not time-reversible. Suppose that it is 
time-reversible. Then, there exists a rotation which pass the variables 
{x,y) to the new variables ("U, f), given by 

(17) M = cos ax — sin ay, t> = sinax + cosay, 

which transforms the axis of symmetry into the line m = 0. In the new 
variables the system becomes ii = P{u,v) and v = Q{u,v). Since this 
system must be invariant by {u,v,t) {u, —v, —t), we must have 

P{u,v) = -P{u, -v), Q{u,v) = Q{u,-v). 

These two equations are satisfied if and only if 

cos a sin a(2 cos a — sin a) = 0, 2 sin^ a — cos^ a = 0. 

Since this system has no solution, the claim is proved. 
Now, consider the following perturbation of system (fT3j) 

HON ^ = i-y + y^){x'^ + y'^-e), 

^^^^ y = (x + 2x2)(x2 + ?/2_^). 

System (fTHj) has also the first integral (fT^ . Consequently, system (fTSj) 
has a center at the origin for all e G M. This center is of linear type 
if £ 7^ 0. Hence, doing the limit e — in system ()18|). we obtain the 
initial system ()15p with a degenerate center at the origin. | 

We have seen in Theorem ^ that the nilpotent centers can be e- 
approximated by systems having linear type centers. Then, applying 
the Poincare-Liapunov method to these linear type centers, and doing 
the limit when e \ we obtain algebraic conditions characterizing the 
existence of a nilpotent center. We shall see this more explicitly in 
Section 0] Now, we shall show that this does not occur for degener- 
ate centers which are not algebraically solvable (proving statement (d) 
of Theorem ^ , and for some classes of degenerate centers which are 
algebraically solvable (proving statement (e) of Theorem El). 

Proof of Theorem\^d): It is known that the problem of determining 
the center conditions at the origin of the system 

X = xp2 — ypi + 4x(x^ + /i?/^)pi, 
y = xpi + yp2 + 4:y{x'^ + fiy'^)pi, 

where pi = x"^ + a^xy + a^y"^ and p2 = ctix^ + a2xy + a^y'^, is not alge- 
braically solvable for some specific values of the parameters, see . If 
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such a system is limit of linear type centers, it would be algebraically 
solvable. So, the proof of statement (d) of Theorem |21 follows. | 

Proof of Theorem\^e): Consider the following cubic homogeneous sys- 
tem 

, . X = P{x,y) = 12Xx^ -9x'^y -20Xxy'^ -25y^ + 9 fiy^, 
^ ^ y = qIx, y) = 9x^ + UXx^y + 25xy'^ - 20A?/^ 

with the monodromy condition that xQ{x,y) —yP{x,y) has no real 
factors. This system has a degenerate center at the origin if and only if 
/i = 0, or A = 0. This follows checking the conditions (i) and (ii) of the 
Appendix which characterize the homogeneous systems having a center 
at the origin. Condition (i) is satisfied by the monodromy condition, 

he) 

and condition fii) is / d9 = 0; i.e. 

-47rA 



v/9/i-25V81/i + 64J-17 - y/MTSTJI J -17 + ^64 + 81/i 



^-17 - ^64 + 81/i(160 - 27/i - 5^64 + 81/i) + 

-17 + ^64 + 81/i(-160 + 27/x - 5y^64 + 81/i)^ = 0. 

It is easy to check that condition (ii) holds if and only if /t = or A = 0. 
Moreover, these centers are algebraically solvable because condition (ii) 
is algebraic. In the case /i = condition (i) is directly satisfied because 
xQ{x,y) — yP{x,y) = {x^ + ?/^)(9x^ + 25?/^). In the case A = 0, the 
homogeneous polynomial xQ{x, y) — yP{x, y) is 9x^ + SAx'^y'^ + 25?/^ — 
9fiy^, and condition (i) is satisfied if and only if /t < 25/9. 
We consider the following perturbation of system (|19|) 



X 



-ey + 12Xx^ - 9x'^y - 20Xxy^ - 25y^ + 9i2y^ + ed, 



^'^^^ y= ex + 9x^ + 12Xx^y + 25xy^ - 20Xy^ + eG2, 

where Gi = Gi{x, y,e), for i = 1, 2, are analytic functions in x, y and e, 
without constants and linear terms in x and y. Applying the Poincare- 
Liapunov method to system (^0]) (see Section El for more details) , we 
obtain that the first Liapunov constant is 

Vi = -8X + eVi, 

where Vi is the first Liapunov constant of the analytic system 

X = -ey + eGi{x,y,e), 
y = ex + eG2ix,y,€). 
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As we must vanish Vi up to first order in e we obtain only the condition 
A = 0. So, the case = and A G M cannot be detected as hmit of 
hnear type centers. | 

Proof of Theorem{^f): We consider the system 

(21) x = -ay^, y = bx^, 

with ab > 0. It is a Hamiltonian system with Hamiltonian 

4 D 

It is easy to check that system (j^U is a (2, 3) -quasi-homogeneous sys- 
tem of weight degree 8 satisfying conditions (i) and (ii) for having a 
degenerate center at the origin, see the Appendix. Another way to see 
that the origin is a degenerate center is noting that the level curves of 
H are ovals. It is easy to check that system (j2T| has only one charac- 
teristic direction, given hj y = 0. 

Now, consider the following perturbation of system (j2H) : 

(22) x = -ay^, y = ex^ + bx^, 

with ae > 0. This system is also Hamiltonian, with 

ay^ ex^ bx^ 
H,{x,y) = ^ + — + —. 

System has a degenerate center at the origin, because the origin is 
surrounded by ovals. This system has no characteristic direction. Now, 
doing the limit e ^ in system ()22|) . we obtain the initial system (j?T|) 
with a degenerate center at the origin and with a characteristic direc- 
tion. I 

The example given in the proof of Theorem Eff ) shows that, in a 
similar way that we can apply the Poincare-Liapunov method to de- 
tect nilpotent centers, in the study of certain degenerate centers with 
characteristic directions we can apply the Bautin method for degener- 
ate centers without characteristic directions (see for instance [31^1) to 
a convenient perturbation of the system with characteristic directions. 

4. The Poincare-Liapunov method for nilpotent and 

degenerate systems 

In this section we illustrate how to apply the Poincare-Liapunov 
method to several families of polynomial differential systems for de- 
tecting nilpotent or degenerate centers. Some of these families have 
been studied recently by other authors with different and more com- 
plicated techniques. First we start studying some nilpotent centers. 
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We note that the simplest nilpotent polynomial centers must be of de- 
gree 3, because there are no nilpotent center for quadratic polynomial 
differential systems, see for instance We consider the system 

(23) X = y + x'^ + k2xy, y = kix^ — x^ . 

We apply to this family the general algorithm, with a general pertur- 
bation and we obtain the following result: 

Proposition 6. System (j23|) has a nilpotent center at the origin if and 
only if ki = k2 = . 

Proof Applying the Poincare-Liapunov method to the perturbed sys- 
tem 

(24) X = y+x'^+k2xy+exGi{x,y), y = —ex+kix'^—x^+exG2ix,y), 
where 

oo oo 

Giix,y)= (^ij^'y\ G2ix,y)= ^ hjx'y^ 

i+j>l «+j>l 

and £ > 0, we obtain the first Liapunov constant 
2 

~(ooi — 3a2o — 2aio6io — ^oi^io ~ ^ii + o,iQk2)£'^ + («02 ~ c^oi'^io)^^ • 

We note that in Vi only appear the linear and quadratic terms of Gi 
and (j2- Vanishing Vi at any order in e we get the necessary center 
condition ki = for system (j^Hj) and for an arbitrary perturbation. We 
obtain also the conditions 6io = ^2/2, aoi = (6020 + 26ii + 6oi^2)/2, 
and ao2 = (6aioa2o — 2aio&ii + aio&oi^2)/2 on the parameters of the 
perturbation. The next Liapunov constant has the form 

V2 = 7 r— ^ ^ [-72k2 + 0(e)] . 

(l + £)(5-2e + 5e2) L ^ 

In the expression of V2 we have contributions of the linear, quadratic, 
cubic, quartic and quintic terms of Gi and G'2. Therefore, the condi- 
tions ki = k2 = are necessary in order that the origin of system (j23j) 
be a center. These conditions are also sufficient as it is explained in 
Remark [T21 | 

We see that it has been sufficient to employ a polynomial pertur- 
bation of degree 5 ir order to determine the necessary and sufficient 
conditions of center for system ()23|). 

Although in Theorem^ the perturbation is unknown, it is surprising 
that with the simple perturbation —ex in y it is possible to obtain 
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the center cases of many families, as it will be shown in the following 
examples. We consider the system 

(25) X = y + Axy + By^, y = —x^ + Kxy + Ly^ . 

Applying the Andreev results |3I we can see that the origin of system 
(|25p is monodromy. 

Proposition 7. System (j^^j) has a nilpotent center at the origin if and 
only ifAB-3L = and AB{A^ - 2K) = 0. 

Proof. Applying the Poincare-Liapunov method to the perturbed sys- 
tem 

(26) X = y + Axy + By^, y = —ex — a;^ + Kxy + Ly^, 
with £ > 0, we obtain the first Liapunov constant 

2e\AB - 3L) 
^ ~ 3 + 2£ + 3£2 ■ 

Vanishing Vi we get the first center condition L = AB/3. Now, we 
compute the second Liapunov constant 

y __ 2e^AB{A'^ - 2K) 
^ ~ ~3(l + e)(5-2£ + 5£2)' 

Vanishing V2 we obtain the second center condition AB{A^ — 2K) = 
0. So, these two conditions are necessary in order that the origin of 
the perturbed system ()26|) be a center. These two conditions are not 
necessary, in principle, for system ()25p. because we must investigate for 
others polynomials perturbations of the form 

X = y+Axy+By'^+exGi{x,y), y = —ex—x^+Kxy+Ly^+exG2{x,y). 

But, in [1^ it is proved that these two conditions are necessary in order 
that the origin of system ()25|) be a center. We remark that in this 
particular system we do not need to take e = in the center conditions 
because they are independent of e. 

Now we prove that these two conditions are sufficient. If A = or 
B = (and consequently L = 0), we have that system (j25|) is reversible 
with respect to {x, y, t) H-i> {—x, y, —t) or {x, y, t) 1— *• {x, —y, —t), respec- 
tively. Therefore, since the origin is monodromy, it is a center. 

If AS 7^ 0, L = AB/3 and - 2K = 0, then the system has the 
analytic first integral 

/ 4 ^ ^ 2 12 12 62 23 , 2 2 „ 3\ 
H = exp{-Ax) [y -JI-J5^- - + + j • 

This first integral can be obtained using the theory of integrability 
of Darboux, see for instance jSOj. In fact, this first integral already 
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appeared in Since the origin is monodromy, by the existence of 
this analytic first integral defined at the origin it follows that the origin 
is a center. We note that in this case the nilpotent center is neither 
time-reversible nor Hamiltonian. | 

We consider the system 

(27) X = —y, y = + ax^ + y{bx^ + cx^). 

Applying Andreev's results 3j we can see that the origin of system (j27j) 
is monodromy. 

Proposition 8. System (P7j) has a nilpotent center at the origin if and 
only if ab = and c = 0. 

Proof. Applying the Poincare-Liapunov method to the perturbed sys- 
tem 

(28) X = —y, y = ex + x^ + ax^ + y{hx^ + ca;"^), 
with £ > 0, we obtain the first Liapunov constant 

v,^ . 

5 + 3£ + 3£2 + 553 

Vanishing Vi we get the first center condition c = 0. Now, we compute 
the second Liapunov constant 

^ _ {2 + 7e)ah 
^~ 128 £2 

Vanishing V2 we obtain the second center condition ah = 0. So, these 
two conditions are necessary in order that the origin of the perturbed 
system (I28|l be a center. These two conditions are not necessary, in 
principle, for system (j77|) because we must investigate, as in the previ- 
ous example, for others polynomials perturbations of the form 

X = —y + exGi{x, y), y = ex + x^ + ax^ + y{bx^ + cx'^) + exG2{x, y). 

But, in |2] it is proved that these two conditions are necessary in order 
that the origin of system be a center. 

Now we prove that these two conditions are sufficient. If a = c = 
or 6 = c = 0, we have that the system is reversible with respect to 
{x,y,t) {—x,y,—t) or {x,y,t) ^— {x, —y, —t), respectively. There- 
fore, since the origin is monodromy, it is a center. | 

Proposition [7| is proved in [1,^ by using Liapunov polar coordinates, 
see jl3j, and computing some generalized Liapunov constants. The 
method developed in [I] is not useful to solve the center problem of 
Proposition see . Proposition |H1 is proved in j2] by using the nor- 
mal form theory and taking into account that a convenient truncated 
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normal form of the nilpotent system is a Lienard system. The method 
developed in this paper solves both problems in a unified form and in 
a more simple way, by computing the Poincare-Liapunov constants of 
a linear center type system. In both proofs we have used the results 
of PP and [2] to prove that the conditions are necessary. This is not 
a restriction of our method because we can apply it with a general 
perturbation. But, in that case it is necessary to make a big amount 
of computations for obtaining the necessary conditions. This is the 
usual amount of computations that appear in the application of the 
Poincare-Liapunov method when the system under study has several 
parameters. 

We consider the system 

(29) X = -y+auxy + ao2y'^ + a3oX^+a2ix'^y+ai2xy'^ + aosy^, y = x^. 

Proposition 9. System (j^Hl) has a nilpotent center at the origin if and 
only if a^o = 0, ao2Ctii + 012 = 0, ci02Oiifl2i = 0, and cio20iiflo3 = 0. 

We consider the system 

(30) X = -y, y = auxy + ao2y^ + a^ox^ + a2ix^y + a^xy^ + ao^y^. 

Proposition 10. System ()3()|1 has a nilpotent center at the origin if 
and only if a2i — 0-02011 = 0, aoa = 0, ao2Ctii'^3o = 0, and 002^11 (3ao2 + 
2au) = 0. 

The proofs of Propositions IHl and ^1 are similar to the proof of Propo- 
sitions [7| and IHl and we omit them. Proposition El is also proved in P] 
and Proposition ITUl is proved in 

As the previous examples show, in some cases, it is sufficient to 
perturb system Q with —ex in y, but there are nilpotent centers which 
are limit of more general perturbations and which cannot be detected 
only with the perturbation —ex in y, as we will see in the following 
example. 

We consider the system 

, . X = P{x, y) = y + xy + {1 — a)y'^ + (1 ~ a)xy'^ — ax^ — ax^, 
y = Qi^-, y) = cy"^ — 2a;^ + cy'^ — 2x^y + (c — 2)x^{l + y). 

Applying the Andreev's Theorem [Hj it is easy to see that system 
has a monodromy singular point at the origin. 

Proposition 11. System (jnU has a nilpotent center at the origin for 
all values of a and c. 
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Proof. System (jHlj) has the following analytic first integral 

which can be determined using the Darboux theory of integrability. 
Therefore, system fl31|) has a center at the origin. | 

Applying the Poincare-Liapunov method to the perturbed system 

y = -ex + Q{x,y). 

with 5 > 0, we obtain the first Liapunov constant 

2A(l + 2a) 
^ ~ 3 + 2e + 3^2 ■ 

Therefore, the first center condition is c(l + 2a) = 0. But, we know 
that system ()31|) has a nilpotent center for all values of a and c. Then, 
it must exist another more general e-perturbation of system (j3ip which 
is a linear type center for all values of a and c. Consider the following 
polynomial perturbed system 

, . X = —ex{ax + ax^)+P{x,y), 

y = -ex{l + (1 - c)x + y + {1 - c)xy) + Q{x,y), 

where P and Q are defined in system (jHH). System has a linear 
type center at the origin because it has the following analytic first 
integral 

H{x, y) = {l + x)-^%l + yy^^ix^ + y^ + ex^). 

Therefore, all the nilpotent centers of system ()3H) are limit of the linear 
type centers of system but not all the nilpotent centers of system 
fl31|) are limit of linear type centers of system ()32|) . This example shows 
that it is not always possible to obtain a nilpotent center as a limit of 
linear type centers with the only perturbation —ex in y, even in the 
case where a local analytic first integral exists. 

Remark 12. Consider the system 

(34) x = y + x'^, ij = — x^ . 

Since this system is time-reversible with respect to the change of vari- 
ables (x, y, t) (—a;, y, —t), and the origin is monodromy (see jHlllj), it 
has a nilpotent center at the origin. But it has neither a local analytic 
first integral, nor a formal first integral defined at the origin, see the 
proof in 110^. 

Consider now the following perturbation of system 

(35) x = y + x'^, y = —ex—x^. 
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As this system is time-reversible with respect to the same change of 
variables, for e > it has also a center at the origin. Therefore, by 
the Linear Type Center Theorem we know that system ()35|) has a local 
analytic first integral H{x, y, e) at the origin. It is possible to compute 
an explicit expression of it given by 



exp 



2 arg(£: + + i{x^ + 2y - e)) (e^ + x^ - 2ey + 2x^y + 2?/^^ 



2^ 

5 



Now, taking the limit when e \ 0, we obtain the first integral 

y) = lim H{x, y, e) = exp 2 arg(x^ + i{x'^ + 2y)) {x^+2x'^y+2y 

e\0 L J 

of system (j34j) . which is not analytic at the origin. □ 

We see in this example that the limit of an analytic first integral 
defined in a neighborhood of the origin can be not analytic. 

In general the study of the nilpotent centers is easier with the algo- 
rithm proposed in this work than applying the results of 0. In our 
case, we have two arbitrary functions Gi and G2, while in the algorithm 
consequence of the results of |8j there are three arbitrary functions, the 
one which appear in the normal form for the nilpotent center and the 
two coming from the change of variables. Moreover, for polynomial 
systems, under the assumptions of Theorem |2l the two arbitrary func- 
tions Gi and G2 of our method are always polynomials and this fact 
does not happen in the algorithm based on the results of [8^ . 

Now we apply the Poincare-Liapunov method to detect degenerate 
centers in a family of polynomial differential systems. 

We consider the polynomial system 

= -a{l + x){x^ - Ay^ -3y^) + iJ,y^, 
^ ^ y = -a{l + y){Ax^ + 3x^ - y^) + Xx\ 

with the monodromy condition a/i > if a 7^ and yuA < if a = 0. 

Proposition 13. System ()36|). with the monodromy condition a/i > 
ifa^O and iiX<Oifa = 0, has degenerate centers at the origin which 
are limit of linear type centers of the form (jT)) with Gi = G2 = if and 
only if fi = X = 0, or a = 0. 



Proof. Applying the Poincare-Liapunov method to the perturbed sys- 
tem 

. . X = ey - a{l + x){x'^ - Ay^ - 3|/^) + 

^ ' y = -ex - a{l + y){Ax^ + 3x'^ -y^) + Xx^, 

with e 7^ 0, we obtain the first Liapunov constant 



e 
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Vanishing Vi we get the first center condition a/i = 0. Now, we compute 
the second Liapunov constant 



Vanishing V2 we obtain the second center condition aX = 0. So, these 
two conditions are necessary in order that the origin of the perturbed 
system fl37|) be a center. Therefore, these two conditions are necessary 
in order that the origin of system be a center which is hmit of 
hnear type centers of the form ((Tj) with Gi = G2 = 0. 

Now we prove that these two conditions are sufficient. If a = 
we have that the system is Hamiltonian and reversible with respect 
to {x,y,t) [x,—y,—t). Therefore, since the origin is monodromy 
(because it has no characteristic directions), it is a center. Now, taking 
the hmit when e ^ 0, we obtain a Hamiltonian system which has a 
degenerate center at the origin. 

If /i = A = it can be shown that system ()36|) has a monodromy 
singular point at the origin. Moreover, it has the analytic first integral 

H{x,y) = {l + x)-\l + yr\x^ + y'), 

defined in a neighborhood of the origin. Therefore, system ()36|) has a 
degenerate center at the origin. We note that this degenerate center is 
neither time-reversible nor Hamiltonian. | 

By the same arguments that for the nilpotent centers of polynomial 
differential systems, the degenerate centers of polynomial differential 
systems which are limit of linear type centers of the form ((Tj) under the 
assumptions of Theorem |21 in fact, are limit of linear type centers of 
the form ((Tj) where the two analytic functions Gi and G2 are always 
polynomials. 

5. On THE CYCLICITY OF NILPONENT AND DEGENERATE CENTERS 

Let p be a center of a polynomial vector field of degree m. The 
cyclicity, Cn{p), of p is the maximum number of limit cycles, taking 
into account their multiplicity, that can bifurcate from the singular 
point p when we perturb it into the class of all polynomial differential 
systems of degree n > m. 

For a linear type center p of a polynomial differential system of degree 
m it is known that if the number of its independent Liapunov constants 
is k, then the cyclicity Cm{p) < k — 1 if the Bautin ideal is radical, see 
for instance jSU]. Moreover, if we perturb a polynomial differential 
linear type center of degree m and cyclicity k — 1 inside the class of all 
polynomial vector fields of degree m, we can get perturbed vector fields 
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with exactly k — 1 hyperbolic limit cycles bifurcating from the center. 
This is due to the relationship between the Liapunov constants and the 
coefficients of the Poincare map near a center. For more details on this 
subject see jHU] . 

As we have seen in the examples, in general, the Liapunov constants 
are not well-defined when £ — >■ 0, see for instance the proof of Propo- 
sitions IHl and El Therefore, we must impose that the limit of the 
Liapunov constants when e ^ be well-defined. If the Liapunov con- 
stants are well-defined when e — >■ 0, then the Poincare map obtained 
by the limit LI = lim^^o gives a formal series which can be not 
convergent at any positive radius. Hence, we must also impose to the 
Poincare map to be convergent in a neighborhood of the origin. Taking 
into account these conditions we can establish the following result. 

Proposition 14. Suppose that the origin of a polynomial differential 
system ^ or 1^ of degree m is a center p, and that this system is 
limit when e ^ of polynomial differential systems of degree n of the 
form or respectively, which have linear type centers p^ at the 
origin. If 

(i) the Liapunov constants of a general perturbation of the same 
degree n of the linear type centers ^ or ^ are well-defined 
when e — >■ 0, and 

(ii) the limit of the Poincare map when e of the general per- 
turbation of the same degree n of the linear type centers or 
((7j) is analytic in a neighborhood of the origin, 

then, 

(a) the cyclicity Cn{p) of the nilpotent center Q is at most 

the cyclicity Cn{pe) of the linear type center (jH)) for all e > 0. 

(b) the cyclicity Cn{p) of the degenerate center is at most the 
cyclicity Cn{pe) of the linear type centers (|7j) for e suffi- 
ciently small. 

Proof. Let (P, Q) be the polynomial vector field of degree m associated 
to the system of the form (0) or Q with a singular point p of center 
type at the origin. By assumptions, the vector field (P, Q) is limit 
when e —>■ of polynomial vector fields (P^, Q^) of degree n associated 
to systems of the form (jH} or ((7|), respectively, which have linear type 
centers Ps at the origin. 

Let {P*, Q*) be a general perturbed polynomial vector field of degree 
n of the vector field (P^, Qs). Taking the limit of {P*, Ql) when e ^ 0, 
we obtain a perturbed polynomial vector field {P*,Q*) of degree at 
most n of the vector field (P, Q). Since the Liapunov constants of the 
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system {P*,Q*) are well-defined when e ^ 0, and the Poincare map 
n = lim^^one of (P*, Q*) is analytic in a neighborhood of the origin, 
we can control the ciclicity of the polynomial vector field {P*,Q*) by 
the Poincare map 11 (with the same restrictions that for linear type 
centers). Moreover, the number of independent Liapunov constants of 
the system (P*, Q*) is at most the number of the independent Liapunov 
constants of the system {P*,Ql). Therefore, the cychcity c„(p) of the 
nilpotent center (j2I) is at most the cyclicity c„(pe) of the hnear type 
center (jH)) for all e > 0, and the cyclicity c„(p) of the degenerate center 
© is at most the cyclicity Cn{pe) of the linear type centers ((Tj) for e 7^ 
sufficiently small. Then the proposition follows. | 

In general the degenerate problems present the more rich structure. 
For instance, when we look for algebraic limit cycles into the qua- 
dratic polynomial vector fields, there is one which is given by a non- 
degenerate algebraic curve (the algebraic limit cycle of degree 2), but 
there are many others (the algebraic limit cycles of degree 4, 5, 6, and 
perhaps others) that are given by degenerate algebraic curves. Here, a 
degenerate algebraic curve is an algebraic curve having singular points. 
For more details about algebraic limit cycles see llj. However, in the 
assumptions of Proposition it is clear that the cyclicity of a non- 
linear type center p which is limit of linear type centers, is not more 
rich than the cyclicity of the linear type centers. Hence, Proposition |3 
follows from Proposition 1141 

6. Appendix: Homogeneous and quasi-homogeneous 

SYSTEMS 

In this appendix we introduce two classes of polynomial vector fields 
having degenerate centers. For more details about them see jT2| and 



We consider polynomial differential systems in of the form 



where P and Q are real polynomials in the variables x and y. We say 
that this system has degree m if m is the maximum of the degrees of 
P and Q. 

If P and Q are coprime homogeneous polynomials of degree m, then 
the centers of systems ()38|) are characterized by: (i) the homogeneous 
polynomial xQ{x,y) — yP{x,y) has no real factors (so m is odd), and 



m 



(38) 



X = P(x, y) 



y = Q{x,y) 



n 
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Here 

f{9) = cos9P{cos6,sm9) + sm9Q{cos9,sm9), 
g{9) = cos 9 Q {cos 9, sin 9) — sin 9 P {cos 9, sin 9). 

Moreover, all the homogeneous centers are global centers; i.e. the pe- 
riodic orbits surrounding the center fulfill all . 

In what follows p and q always will denote positive integers. 

We say that the function H{x,y) is {p,q) -quasi-homogeneous of 
weight degree m > if H{Fx, I'^y) = V^H{x, y) for all / G M. 

We say that system (jHH|) is (p, q) -quasi-homogeneous of weight de- 
gree m > if P and Q are (p, -quasi-homogeneous functions of 
weight degrees p — 1 + m and q — 1 + m, respectively. Note that the 
(1, l)-quasi-homogeneous systems of weight degree m are the classical 
homogeneous polynomial differential systems of degree m. We note 
that if system (jHHj) is {p, g) -quasi-homogeneous, then the differential 
equation dy/dx = Q/P (another way to write system (jHHj) ) is invariant 
by the change of variables {x,y) —>■ {PxJ'^y). 

If P and Q are coprime, then the centers of the {p, g) -quasi-homogeneous 
systems of degree m are characterized by: (i) the {p, g)-quasi- 
homogeneous polynomial pxQ{x,y) — qyP{x,y) has no real factors, 
and (ii) 

r ^^^^d9-Q 

Jo W) 

Here 

F{9) = Cs^''-^9P{Cs9,Sn9) + Sn^P-^9Q{Cs9,Sn9), 
G{9) = pCs9Q{Cs9,Sn9) -qSn9P{Cs9,Sn9), 

and Cs 9 and Sn 9 are the (g, p)-trigonometric functions. Moreover, all 
the {p, g)-quasi-homogeneous centers are global centers. 

We recall that the {p,q) -trigonometric functions z{9) = Cs9 and 
w{9) = Sn^ are the solution of the following initial value problem 

i=-ti;2p-\ ^ = ^2g-l^ ^(o)=p"^, w(0) = 0. 

It easy to check that the functions Cs 9 and Sn 9 satisfy the equality 

pCs^'^^ + gSn^*'^ = 1. 

For p = g = 1 we have that Cs = cos 9 and Sn ^ = sin 9\ i.e. the (1, 1)- 
trigonometric functions are the classical ones. The functions Cs6' and 
Sn^ are r-periodic functions with 
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